Draft version February 5, 2008 

Preprint typeset using WT^^ style cmulateapj v. 6/22/04 



ACCRETION-EJECTION INSTABILITY, MHD ROSSBY WAVE INSTABILITY, DISKOSEISMOLOGY, AND 

THE HIGH-FREQUENCY QPO OF MICROQUASARS 

Michel Tagger 

Service d'Astrophysiquc, (UMR AstroParticules et Cosmologic), CEA Saclay 
91191 Gif-sur-Yvette, France 



o 
o 

(N 

< 



> 

m 
(N 

oo 
O 

o 

:^ 

o 



X 



Peggy Varniere^ 

LAOG, Universite J. Fourier (UMR5571) France 
Draft version February 5, 2008 

ABSTRACT 

We present a possible explanation for the high-frequency Quasi-Periodic Oscillations of microquasars 
by an MHD instability that combines the physics developed, in different contexts, for the Accretion- 
Ejection Instability, the Rossby-Wave Instability, and the normal modes of diskoseismic models (which 
rely on the properties of the relativistic rotation curve in the vicinity of the Marginally Stable Orbit). 
This instability can appear as modes of azimuthal wavcnumbers m = 2,3, . . .that have very similar 
pattern speeds w/m, while the m = \ mode, which would appear as the fundamental of this discrete 
spectrum, is less unstable. This would readily explain the 2:3 (and sometimes higher) frequency ratio 
observed between these QPO. 

These instabilites form eigenmodes, i.e. standing wave patterns at a constant frequency in the disk; 
they arc strongly unstable, and thus do not need an external excitation mechanism to reach high 
amplitudes. Furthermore, they have the property that a fraction of the accretion energy can be 
emitted toward the corona: this would explain that these QPO are seen in a spectral state where 
Comptonized emission from the corona is always present. Their existence depends critically on the 
existence of a magnetic structure, formed by poloidal flux advected in the accretion process, in the 
central region between the disk and the black hole. 
Subject headings: accretion — black hole physics-X-rays: 
MHD— Instabilities 
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1. INTRODUCTION 

Although many models have attempted to explain 
the Quasi-Periodic Oscillations (QPO) of X-ray bina- 
ries, none of them has gained yet a wide acceptance, 
and many have been found in contradiction with the 
observed frequencies and properties of the QPOs. The 
latest addition to this mystery has emerged in the last 
few years: it is the observation of high-frequency QPO 
(HF-QPO) in microquasars, where they are observed 
at different fr equencies in a 2:3 (and sometim es higher) 
ratio (see e.g. iMcClintock fc R,emillardll200fitl . 
On the other ha nd, we have in rece nt years 
lIRodnguez et a,1J lVa,rniere et ahl 

iTagger et al.l I2004D proposed to explain the low- 
frequency QPO (LF-QPO), observed when the sources 
are in different spectral states, by an MHD instability 
of magnetized accr etion disks, the Accret ion-Ejection 
Instabihty (AEI, iTagger fc Pellaa EMI- Besides 
producing a frequency in the right range, and other 
properties which support this explanation, it has the 
advantage of being a normal mode, i.e. to produce 
long-lived patterns at a quasi-stationary frequency, and 
of being an instability, i.e. to grow spontaneously to 
high amplitudes if the disk is in a favorable magnetic 
configuration. The AEI is formed of a combination of 
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spiral and Rossby waves propagating in the disk, and 
coupled by differential rotation. Its name relates to the 
fact that, if the disk is covered by a low-density corona, 
a sizable fraction of the accretion energy extracted from 
the disk can end up in Alfven waves emitted by the 
inst ability in the co rona, where they might energize a 
jet ijVarniere fc Tagg er 2002). This might correspond to 
the fact that the LF-QPO is seen when the source is 
in the 'low-hard' state, where the disk emission is very 
faint and dominated by a strong and hard emission from 
the corona. In this sense, although it exists already in 
an infinitely thin disk in vacuum, and we are able only 
to treat a low-density corona as a perturbation, the AEI 
must be viewed as a three-dimensional instability of the 
disk-corona system. 

In this paper we show that the HF-QPO, observed 
in a few microquasars when they are in a Very High 
(or steep-power-low) state, may be explained by an 
instability, due to the specific conditions near the 
Marginally Stable Orbit (MSO) of a black hole disk, 
and which can be described as a hybrid of the AEI 
and of a different, but closely related instability: an 
MHD form of the Ros s by-W ave Insta bility (RWI: 
ILovelace fc Hohlfeldl ITOTSt iLcwela cc ct aj] 119991) . Or, 
to present this differentl y, we show that the mod es of 
diskoseismologic models ()Nowak fc Wagoneil I1991D are 
strongly unstable by the AEI-RWI mechanism. All 
the modes with azimuthal wavenumber m > 2 have 
very close pattern speeds w/m and growth rates so 
that, depending on the conditions, any of them can be 
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expected to dominate, readily producing frequencies in 
the ratio 2:3:4... On the other hand the m = 1 mode 
(which would appear as the fundamental in this discrete 
spectrum) is always much less unstable. Although the 
instability exists already in unmagnetized disks, it is 
much stronger when the disk is thread ed by a vertical 
magn etic field. We have recently shown ijTaeeer fc Melial 
|2006|) that the RWI could explain the quasi-periodic 
behaviour observed during the flares of the central black 
hole of the Galaxy in Sgr A*, when a blob of gas from 
the environment is captured in the disk at a few to a 
few tens of Schwarzschild radii. Here on the other hand 
the instability is due essentially to the properties of the 
relativistic rotation curve, in the vicinity of the MSO. 
Our paper is organized as follows: in section |21 we will 
present the basic physics of Rossby waves, and the 
conditions that can make them unstable (becoming 
the RWI) in a differentially rotating disk. We will also 
present the connection between the RWI and the AEI. 
which share much of their physics. In section 12.21 we 
will discuss the properties of diskoseismologic modes 
that exist in the specific conditions near the MSO of 
relativistic disks, and show that they are unstable by 
our mechanism. Section O will be devoted to numerical 
results, first via a linearized analyis and then by non- 
linear MHD simulations which lead us to discuss the 
magnetic conditions near the inner radius of the disk. 
The last section will be devoted to a discussion of these 
results and their astrophysical implications. 

2. MAGNETIC ROSSBY WAVE AND INSTABILITY 

2.1. The 'classical' RWI 

The RWI has been found and discussed in differ- 
ent contexts of differentiall y rotating disks. In g alacti c 
disks, it was introduced bv lLovelace fc Hohlfeldl l)1978D . 
who called it Negative-Mass Instability. They showed 
that a disk presenting an extremum of a certain quan- 
tity C (later dubbed vortensity, see below for its defini- 
tion) was subject to a local instability of Rossby vor- 
tices. The requirement of an extremum is similar to 
that giving rise to th e Kelvin-Helmholtz ( KH) instabil- 
ity of sheared flows (jPrazin fc R,eidl [20041) . so we may 
view this as an application of the KH inst ability to disks 
in diff erential rotation; more recen tly, iLovelace et al.l 
l|1999t> (see also lLi et al]l200nl l200l developed the the- 
ory and numerical simulation for accretion disks and 
renamed it the Rossby Wave Instability (RWI). In the 
meantime it had been re-discovered in numerical simu- 
lations of galactic disks, and the theory re derived, by 
iSellwood fc Kahl l(T99l|) : iDubruUe fc Zahl |T99l|) also 
discussed it (as a non-linear instability) in t he context 
of accretion disks. iVarniere fc Taggejl (|200fitl used it to 
discuss a possible scenario for accretion and planet for- 
mation in the 'Dead Z one' of protostellar disks. Finally 
iTagger fc Melial l)2006|) pointed to the MHD form of the 
instability, and showed that it might explain the quasi- 
periodic modulation seen during the flares in Sgr A*. 
In isothermal, unmagnetized disks, the quantity C is the 
specific vorticity averaged across the disk thickness, 

W 

>C - ^ (1) 



where 



and S is the disk's surface density, its rotation fre- 
quency, and K the epicyclic frequency, defined by 

= 4n^ + 2nn'r^—^(nr^) . (3) 

r dr 

The prime notes the radial derivative, and one notes that 
K = O in keplerian disks. In ma gnetized disks C is re- 
placed by (|Tagger fc PellatilT999|) 
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Although the RWI has been already discussed in the 
above-mentioned references, we present here its basic 
physics for the sake of completeness. It results from 
two essential properties of waves in differentially rotating 
disks: 

• Rossby waves are most often studied in fluids pre- 
senting a solid rotation, but a gradient of vorticity, 
such as planetary atmospheres: the rotation is that 
of the planet, while vorticity has a gradient with 
latitude. They propagate perturbations of vortic- 
ity in the flow. Their dispersion relation is 
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where x is the direction of the gradient (latitude, 
or radial in what follows), y its normal (longitude, 
or azimuthal) and W is the latitudinal gradient of 
vorticity. 

This dispersion relation shows the important prop- 
erty that Rossby waves can propagate only in one 
direction, or equivalently produce vortices of only 
one orientation, cyclonic or anticyclonic, depend- 
ing on the sign of W' . Differential rotation results 
in a D oppler shift of the dispersion relation (jTaggeil 
I2OOI . giving in a disk: 
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where m is the azimuthal wavenumber. This ap- 
plies in unmagnetized, isothermal disks of constant 
surface density (so that W appears here rather 
than C) whereas, as shown in the above-mentioned 
works, the vorticity gradient W that appears in the 
numerator must include additional terms in more 
complex conditions (gradients of density, giving £, 
and magnetic field strength, giving £g, polytropic 
equation of state, etc.). 

• The fact that depends on r, while as discussed 
above the right-hand side of equation has a def- 
inite sign (choosing m positive for dcfiniteness) 
means that, depending on the sign of W' (or the 
gradient of C or Cb depending on the physics de- 
scribed) , a Rossby wave can propagate only on one 
side of its corotation radius, defined as the radius in 
the disk where the left-hand side of equationElvan- 
ishes. Since in general the energy of waves is posi- 
tive beyond corotation, and negative inside it, this 
also means that the wave can have only either pos- 
itive or negative energy, contrary to density waves 
that can propagate on both sides of corotation and 
have accordingly positive or negative energy. 
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• The dispersion relation also shows that Rossby 
waves in disks propagate in an annulus starting 
at the corotation radius, contrary to spiral den- 
sity waves which have a forbidden band centered at 
corotation but can propagate on bot h sides beyond 
this band. It was however shown ()Taggeill2001() 
that a Rossby wave is usually sheared away by dif- 
ferential rotation before it has completed even one 
oscillation cycle, unless some mechanism (such as 
the instability mechanisms discussed here) regen- 
erates it continuously. 

Most disk instability mechanisms rely on an exchange 
of energy (and the associated angular momentum) be- 
tween waves of positive and negative energy: this is 
the case, for instance, of the Swing amplification of spi- 
ral waves in galaxies (where waves are driven by self- 
gravity), w hich also applies to the P apaloizou-Pringle 
instability ijPapaloizou fc Pringldll985|l (where they are 
driven only by the gas press ure), and also of th e MHD 
instability first described bv iTagger et all l)199(]() where 
they are driven by the Lorentz force. In these three cases 
a spiral wave propagating inside corotation can be am- 
plified by coupling to another spiral propagating beyond 
corotation, with which it exchanges energy and angu- 
lar momentum. Amplification results from the fact that 
the wave inside corotation can increase its amplitude 
by transfering energy and momentum to the wave be- 
yond corotation. The spirals can also couple to Rossby 
waves near corotation. This process, known as corota- 
tion resonance, is stabilizing or destabilizing depending 
on the gradient of £ which decides whether the Rossby 
wave has a positive or negative energy. Given the form 
of £, this is usuall y stabilizing for spiral modes in un- 
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magnetized 

[197 6: Papaloizqu fc Pringlc JL9871; Nar avan et al.l 



Mark! 
19871 



Papaloizou fc Linlll989(l . unless the radial density gradi- 
ent is unusually strong. On the other hand, since the gra- 
dients act differently on £3, it has been found to give rise 
to the Accretion-Ejection Instability in disks threaded by 
a vertical magnetic field ( Tagg er fc Pclla t 199^. 
Modes amplified by the Swing or by the corotation res- 
onance depend on the existence of an inner disk edge 
where the spiral waves can be reflected. This allows them 
to travel back and forth between the inner edge and the 
corotation radius, forming standing patterns (i.e. normal 
modes) when their frequency is such that they verify an 
integral phase condition. 

The situation changes however when C (or Cjs) has an ex- 
tremum in the disk. Let us first consider the case where 
it has a minimum: then according to the dispersion re- 
lation © a Rossby wave whose frequency lies precisely 
at this minimum can propagate both beyond it, where 
C is positive, and within it, where it is negative. The 
above-mentioned works have shown that this allows nor- 
mal modes, i.e. standing wave patterns trapped near the 
minimum, and that these modes are unstable, i.e. grow 
exponentially by exchanging energy between their pos- 
itive and negative energy components. When C has a 
maximum, one also finds similar unstable modes, with 
now a corotation radius offset from the maximum of C: 
these modes are localized where Rossby waves can prop- 
agate, i.e. between the corotation radius and the maxi- 
mum of C. This instability mechanism was used recently 
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Fig. 1. — Rotation frequency Q, (solid), epicyclic frequeney k 
(dots) and W = K^/2f2 (dashed) in the Paczinsky-Witta pseudo- 
Newtonian potential. 



bv iTagger fc Melia[ l[2006j) to explain the quasi-periodic 
behavior observed during the flares of Sgr A*, by the 
capture in the disk of a fragment of stellar wind which 
creates at a few to a few tens of gravitational radii an 
extremum of density and magnetic field, and thus of C13. 
The RWI thus establishes a standing vortex in the coro- 
tation region. Since it essentially exchanges gas (and 
magnetic flux) from inside corotation with gas from out- 
side, it might be viewed as a form of the classical MHD 
interchange instability. However here the main energy 
source driving the instability is differential rotation (i.e. 
ultimately gravitation) , much stronger than the magnetic 
fleld gradient responsible for the interchange mode, un- 
less the disk becomes magnetically (rather than ccntrifu- 
gally) supported against gravity. 

2.2. Diskoseismologic modes 

In all of these earlier applications, the extremum of 
C was mostly considered as due to an extremum in the 
radial density profile. However, starting from a very 
different point of view — the seismol o gy of relativistic 
accretion disks — iNowak fc Wagoiiei[ (|199H) described 
modes associated with the same extremum of £, but due 
this time to the behavior of k: specifically, the fact that 
the relativistic rotation curve in the disk of a compact 
object is characterized by the existence of a Marginally 
Stable Orbit (MSO), where k vanishes. Inside this 
orbit (at r^so = QGM/c^ for a Schwarzschild black 
hole), K is imaginary, meaning that the motion of fluid 
particles is unstable and they rapidly spiral toward the 
compact object. Moving toward larger radii beyond the 
MSO K first increases, but then decreases again since 
K w fi ~ r""^/^ in the Keplerian portion of the disk. 
Thus K and W have maxima, and usually so does £, as 
shown in figure ^ slightly beyond t he MSO. 

The non-axisymmetric g-modes of INowak fc Wagoneil 
which are trapped in this extremum, are thus similar 
to the RWI, but it was not realized at that time that 
this similarity made them unstable (and very strongly. 
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as will be seen below) by the same mechanism found 
in the above-mention ed works. In wh at follows we will 
use the definition of iLi et all l)2003j) . who classify as 
g-modes the ones for which 'most of the action takes 
place near corotation', i.e. the RWI (whereas the AEI 
would be classified as a p-modc). In that work it was 
shown that the non-axisymmetric g modes which have 
nodes in their vertical structure are always stable. We 
are interested here in the fundamental g-modc, the 
one without nodes, i.e. the one whose structure is 
essentially constant across the vertical extension of the 
disk. This allows us to consid er an infinitely thin disk 
in vacuum ijTagger et al.lll992l). as was the case e.g. for 
galactic spiral m o des (iGoldreich fc Lvnden-Belll 119651 
iSellwood fc KahnI Il991|) before 3D simulations were 
possible. On the other hand, this approach forbids us 
to consider the Magne to-Rotational Instability (MRI, 
iBalbus fc Hawlevlll99H) . which depends on variations in 
the vertical direction: technically, when one studies the 
vertical structure of the modes in a disk with a realistic 
density stratification, one finds a discrete set of modes 
with Uz = 0, 1, 2, 3... nodes (this can be replaced by 
a vertical wavenumber when becom es large. The 
mode s we consider correspond to = (|Tagger et al.l 
Il992j) . whereas the MRI corresponds to riz > 1. Only 
fully three-dimensional simulations with the proper 
magnetic topology would permit to accomodate both 
types of instabilities. 

3. NUMERICAL RESULTS 

In the present work we are mostly interested in the 
MHD form of the RWI, in a disk threaded by a vertical 
magnetic field, for a number of reasons: one is that accre- 
tion disks are magnetized, and that accretion from larger 
radii to the inner regi on should rely on the M agncto- 
Rotational Instability (|Balbus fc Hawlevlll991l ^. A sec- 
ond reason is that in general we find that magnetic 
stresses make the modes significantly more unstable than 
in the unmagnetized case. 

A last but key reason is that in previous work 
(jVarniere fc Taggeill2002|) dedicated to the AEI we have 
derived a very important property of Rossby waves in 
that configuration: they represent a vortical motion ap- 
plied to the fiuid, and thus also to the footpoints of the 
magnetic field lines threading the disk. In an unmagne- 
tized disk the linear physics would stop here, and one 
would have to turn to non-linearities to consider the ul- 
timate evolution of the energy stored in the Rossby vor- 
tices, presumably when the perturbed motions become 
comparable with the sound speed (|Naravan ct al. 1987). 
On the other hand, if the disk is threaded by a vertical 
magnetic field, this vortical motion will propagate along 
the field lines into the corona of the disk as Alfven (or 
more precisely Rossby- Alfven) waves. Our result showed 
that these waves could carry to the corona a significant 
fraction of the accretion energy extracted from the disk. 
This provides naturally a channel to energize the corona 
and maybe form a jet or a wind. We will return to this 
point in our final discussion. 

Since the physics and structure of the modes is essentially 
similar to the unmagnetized case we will not repeat here 
the linear theory. We will present in the next section 
numerical results, first solving the linearized set of equa- 




FlG. 2. — Pattern speed Lj/m (solid) and growth rate 7 (dashed) 
of the Hnear eigenmodes. Top, as a function of /3 = Sirp/B'^; bot- 
tom, as a function of m for (3 = 5. The results are normalized to 
the rotation frequency at the MSO. 



tions for the unstable modes, and then turning to fully 
non-linear MHD disk simulations. 

3.1. Linear analysis 

We first use the same procedure as in iTagger fc PelO 
(|1999|) . to which we refer for a more detailed description: 
we write the linearized MHD equations for perturbations 
in the disk, and project them on a radial grid of points 
equally spaced in s = Inr. The magnetic field in the vac- 
uum surrounding the disk is described by a magnetic 
potential, whose source is the currents in the disk. The 
equilibrium magnetic field is taken normal to the disk. 
For simplicity, since the gas pressure plays little role, we 
use an isothermal equation of state, and at the MSO we 
use a reflecting boundary condition, u.,. = 0. We defer 
a discussion of this boundary condition to the next sec- 
tion, where this will be compared with other boundary 
conditions. At the outer radial boundary we use a radi- 
ation condition (i. e. waves are allowed to freely escape), 
implemented as in ITagger fc Pella^ ljl999() by solving on 
a complex s axis. This also improves the numerical con- 
vergence in the corotation region. 

The only difference is thus that here we use the pseudo- 
newtonian potential of Paczinsky-Witta: 



This mimicks the relativistic rotation curve by allow- 
ing for the existence of an MSO, where k vanishes, at 
I'M so — 3^5 = QGM/(? for a Schwarzschild black hole. 
Figure ^ shows the radial profiles of and k in this 
potential. The resulting discretized equations form a 
(4nr X Aur) matrix, whose eigenvectors and eigenvalues 
give respectively the spatial structure of the perturbed 
quantities and the corresponding frequencies. Physically, 
these eigenmodes are standing wave patterns that have 
an exponential variation with time, defined by a real fre- 
quency (giving an azimuthal pattern speed Lo/m) and a 
growth rate. Among the 4nr eigenvalues one easily sepa- 
rates by various means (essentially a convergence study) 
those of numerical origin, that result from the discretiza- 
tion, from the physical spectrum, in which we consider 
only the fundamental, most unstable mode. 
For the results presented here we use simple power- 
law profiles for the surface density and magnetic field, 
Y, r~^l'^ and B ^ 7'"^, and a constant temperature. 
The only free parameters are then the sound velocity Cg 
and the magnetization parameter /? = %Trp/B^. We take 
varying values of /3, and Cg ~ .05 rQ{rMso)j permitting 
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a clear separation of scales with the rotation velocity. 
The behaviors we obtain are very weakly sensitive to the 
value of cs, since thermal pressure plays only a minor 
role in the physics of the modes we consider (on the other 
hand decreasing cs would decrease the wavelength of the 
most unstable mode, in the unmagnetized case). Since 
the physics is very localized near the extremum of Cb, 
different profiles make very little practical difference. 
Figure [3 shows the frequency and growth rates of modes 
as a function of m for (3 = b, and as a function of [3 for 
m = 1 to 4. These results show that 

• as expected, the modes have very close azimuthal 
pattern speeds {oj/m), that correspond to corota- 
tion radii very near the inner radius. They will 
thus appear very near multiples of a fundamental 
frequency. One has to be cautious about its precise 
value, since it depends on our assumption (mod- 
eling the relativistic rotation curve by a pseudo- 
ncwtonian potential) and also on the spin of the 
black hole. 

• although the instability already exists when /? be- 
comes very large (in the limit /? — > oo one gets 
the hydrodynamical case of Lovelace and cowork- 
ers), one gets much higher growth rates ast f3 ap- 
proaches 1 (equipartition between the gas pressure 
and magnetic field pressure). 

• the growth rate of the m ~ 1 mode is always 
smaller than for m > 1. Modes with m > 1. on 
the other hand, have very similar growth rates, so 
that we can expect minute differences in the disk 
physics, not included here, or in non-linear effects 
to select any value of m > 2 as the mode that 
will dominate oscillations of the disk. This readily 
makes this an excellent candidate to generate the 
observed HF-QPO at frequencies 2, 3, and some- 
times higher multiples of an unobserved fundamen- 
tal frequency. 

• The growth rates we find are very strong, giving 
exponentiation time of the order of the rotation 
period. This is due to the very strong gradient of 
Cb near the inner edge. 

In the next section, we present the results of MHD simu- 
lations that will help us explore further these questions. 

3.2. Numerical simulations 

We perform numerical simulations with the code pre- 
sented in^aunt & Taeeer (2001), to which we refer for 
a detailed description. As for the linear computations of 
section [3.11 it studies the evolution of an infinitely thin 
disk in vacuum, and the magnetic field outside the disk is 
derived from a magnetic potential related to the vertical 
component of the field at the disk by a Poisson equation. 
This makes the problem very similar to simulations of 
flat, self-gravitating disks and allows us to use me thods 
(the Poisson kernel of Binn^y ^ Tremaiii3 l)1987l) . and 
the FARGO algorithm of lMassetl ^200($\ ) developed in 
that context. The code uses a 2D cylindrical (r, 0) 
grid, logarithmically spaced in radius to allow a good 
resolution in the inner disk region, while the outer 
boundary is kept far enough to avoid the influence of 
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Fig. 3. — Azimuthally averaged profiles of surface density, radial 
velocity, magnetic field, 13, C and Cs in our first simulation, with 
a free flow boundary condition at -Gtmso- Radii are scaled to 
f'MSO- Only the inner region is shown, although the simulation 
extends to 15rj\,/so. 



waves reflect ed there. A cons ervative scheme of the type 
described bv lStone & Norman. (4992; ) is used to advance 
the fluid properties in time. The resolution is typically 
taken to be 256x128 cells in r and although higher 
resolution tests have been made for convergence s t udies . 
The main changes to the code of lCaunt fc Taggeil (|2001fl 
are the use of the pseudo-newtonian potential discussed 
above, and of various boundary conditions at the inner 
edge of the disk, as presented in the following sections. 

3.3. Free flow boundary condition 

We first present simulations that allow the gas, once 
it has passed the MSO, to flow freely toward the black 
hole: it thus forms a 'plunging region', as seen in 3D sim- 
ulatio ns that use the same hyp othesis ijHawlev fc Kroli^ 
l200lt iMachida fc Matsumotd 1200^ . This is imple- 
mented by fixing the inner grid point at typically r = 
■ 6 rMSOi with 30 grid points between this and the MSO, 
and allowing the gas to flow inward at this boundary. A 
plunging region is established self-consistcntly in ~ one 
rotation time, and doesn't vary much afterward. Typi- 
cally a sonic point is passed smoothly inward from the 
MSO, and the radial velocity at the inner point is ~ Scg. 
This also creates an outgoing pulse of density propagat- 
ing outward at the sound speed, so that a few more rota- 
tion times are needed to establish quasi-stationary pro- 
files in the inner few MSO radii. The resulting radial 
profiles are shown in figure |21 

We might use these profiles as input in the numer- 
ical procedure of section 13.11 for a linear solution, but 
find it unnecessary since at very low amplitudes the code 
provides simultaneously and reliably a linear solution for 
modes at different m values. Also, the maximum of C 
and Cjs near the MSO keeps evolving as the inner disk 
edge is slowly eroded by accretion caused by the insta- 
bility. We could use viscosity (representing small-scale 
turbulence, presumably due to the MRI) to provide a 
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Fig. 4. — Radial velocity in our first simulation with free-flow 
boundary condition, showing the strong m = 1 mode appearing as 
a one-armed spiral, for enhanced visibility the velocity is plotted 
both as a fake third dimension and as colors, and the disk is seen 
at an oblique angle. 




Fig. 5. — Surface density in the disk. 



continuous inflow that would stop this erosion, but this 
would require additional modeling which is unnecessary 
at this stage. In practice, we stop the simulation when 
non-linearities become important, since they would prob- 
ably be strongly affected by our thin-disk assumption and 
by the limited physics in the code, and at the present 
stage we restrict our discussion to the linearized physics. 

In these conditions, we do not expect modes, such 
as the AEI (or the Papaloizou-Pringle instability in the 
unm agnctized c ase) that depend on a reflective bound- 
ary (iBlaesHigSTj) . On the other hand, as shown in figure 
El the radial profile of k creates a sharp maximum of 
Cb just beyond the MSO, so that we can expect (and do 
obtain) a relativistic g-mode, unstable by the RWI mech- 
anism. For all 'reasonable' radial profiles of density and 
magnetic fiux, and values of the physical parameters, we 
find that the m = 1 mode is always dominant. Figures 01 
and [S] show surface plots of the resulting radial velocity 
and surface density, showing the strong 1-armed spiral. 
Detailed analysis shows that this mode has a corotation 
radius at « 1.3 tmsOi is trapped in the extremum of 
Ci3^ . It is thus the expected non-axisymmetric mode of 

^ This is exactly where we expect it given the profile of Cb- 
in the WKB limit, a Rossby wave can propagate inside corotation 
when the gradient of Cb is negative. Thus here it can propagate (in 
both radial directions, permitting the formation of a standing mode 
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diskoscismology, unstable by the RWI mechanism. But 
simulations exploring a wide range of parameters and ra- 
dial equilibrium profiles always show a strongly dominant 
m = 1 mode, ruling out this mechanism as an explana- 
tion for the observed HF-QPO. 

3.4. Reflecting boundary condition 

On the other hand, this boundary condition that lets 
the gas and its magnetic flux flow freely to the black 
hole may not be the relevant one. This question actu- 
ally reduces to whether or not the disk carries a non- 
vanishing vertical (poloidal) magnetic flux, and what 
happens to that flux when it has b een advected to the 
inner disk edge. It has been argued ijLubow et alJll99^ 
that turbulent resistivity might prevent the poloidal fltix 
to be dragged with the gas in the disk. This relies how- 
ever on rather crude assumptions and the result can- 
not be taken as granted. For instance, axisymmet- 
ric si mulations of magnetized accretion-ejection struc- 
tures (jCasse fc KeT)Den^l2002() . performed with hypoth- 
esis which should be favorable to this effect (turbulent 
magnetic diffusivity much stronger than viscosity), do 
show the poloidal flux advected with the gas. This is 
obtained in axisymmetric simulations, and thus does not 
even rely on the non-axisymmetric me chanism recently 
discussed bv lSpruit fc Uzdenskvl ()2005j) . 
Most current three-dimensional MHD simulations of disk 
turbulence, in newtonian or gencral-relativistic gravity, 
consider disks with no net poloidal flux. This avoids a 
number of numerical and physical problems associated 
in particular with the boundary conditions, and with the 
initialization of the code. But these simulations have 
thus far found no QPO and, although they find jets from 
some form of the Blandford & Ziiaick (19X3) mechanism, 
these jets are not colliniated. 

Most MHD models of colliniated jets powered by the 
disk (rather than the black hole spin) do require a net 
poloidal flux both to accelerate the gas and to coUimate 
it. This is the case for centrifugally o r magnetically 
drive i i iets (iBlandford fc Pavnd 119821: iLovelace et alJ 
119871 iPelletier fc Pudritzlll992K It is also required for 
the AEI, which provides a possible explanation for the 
Low-Frequency QPO in that configuration. Since these 
models require a magnetic field of the order of equipar- 
tition (/3 ~ 1), it is likely that the poloidal fiux that 
is established, in present 3D simulations, near the inner 
edge of the disk is too weak (/3 ~ 10) and localized to 
provide coUimation. 

If poloidal magnetic flux is advected with the gas, one has 
to consider that it must pile up in the central cavity be- 
tween the disk and the event horizon: this leads us to dis- 
cuss the black hole magnetospheric structure that holds 
this flux, involved i n the General-R elati vistic MH D pro- 
cesses of Blandford fc Znaiekl ()1977f) and lPenros'3 l|r969). 
Although thr ee-dimensional n umerical models are avail- 
able (sec e.g. iDe Villiers et al.. .200511) . they are still ob- 
tained in a configuration where the disk has no net 
magnetic flux. The most advanced axisymmetric sim- 

pattern) between the maximum of Cb at r=l.l and the corotation 
radius at 1.3. If Cb had a minimum, the corotation would be right 
at the minimum, and Rossby waves could propagate on both sides 
of it be cause of the inverse gra dients of Cg- This is for instance the 
case m ITagger fc Melial S20M\ . where a minimum of Cb is created 
by a maximum of density. 
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Fig, 6, — Azimuthally averaged profiles of surface density, radial 
velocity, magnetic field, /3, C and Cg in our second simulation, 
with a reflecting boundary condition at the MSO. Radii arc scaled 
to rjvfSO- 



ulation ()Komissarovll2005|) . involving a Kerr Black Hole 
embedded in an external magnetic field, does show the 
Blandford-Znajek process at work, though apparently 
not efficiently enough to produce a relativistic jet. 
This central magnetospheric structure can only be due to 
currents outside the Black Hole; it is commonly consid- 
ered that they form a force-free structure in the magne- 
tosphere, and a current ring which must sit at the inner 
edge of the disk since it must be tied to some mass in 
order to avoid rapid expansion. This can be seen for in- 
stance in Q<!omissarov (2005) where a quasi-steady state 
is obtained after a few rotation times, and where the sub- 
sequent magnetic flux evolution is discussed in terms of 
slow resistive processes. If this is the case one can expect 
that, on the fast time scale of the modes we consider, this 
magnetospheric structure will oppose the penetration of 
the gas from the disk and thus effectively act as a reflect- 
ing boundary: it would oppose the fast flow associated 
with a plunging region, discussed in section rOl or with 
the waves we consider here; but it would not preclude a 
gas flow on the slow time scales of reconnection and/or 
accretion. 

A complete solution of this problem is beyond the pos- 
sibilities of present simulations, and we simply replace 
here this central magnetospheric structure by the same 
boundary condition (vanishing radial velocity at the in- 
ner radius) as in the linearized computation of section 

EH 

The most natural hypothesis on the location of this 
boundary (the disk/magnetosphere transition) is the 
MSO: both because this is where the gas dynamics any- 
how undergoes a transition (more precisely, the inner 
disk radius cannot be smaller than the MSO without 
a very strong magnetic support, of the order of grav- 
ity rather than thermal pressure as assumed here), and 
because of observational results: in the high states of 
microquasars, spectral fits indicate that the color radius 
consistently reaches a minimum and steady value. The 
relation between this and the actual inner radius of the 




Fig. 7. — Radial velocity in our second simulation, showing a 
dominant m = 3 mode. 




Fig. 8. — Surface density in the disk. 



disk must be taken carefully, since it is model dependent; 
but the evidence of this minimum value certainly tells us 
that the inner radius also reaches a minimal value, which 
in all likelihood must be the MSO. 

The results of these simulations globally confirm the re- 
sults of section 13,11 starting from very low amplitude 
noise (from which the unstable modes grow), we can 
accurately measure the growth rates of the different m 
modes, and find them in good agreement with the lin- 
ear results, both for the radial structure of the perturba- 
tions and for the frequencies and growth rates. Changing 
the initial noise component, or adding some viscosity to 
mimic the action of the MRI (supposed to act at smaller 
scale if we stay above its threshold /3 = 1) allows the 
dominant mode to vary but the difference in growth rate 
with the m = 1 mode is strong enough that only artificial 
initial conditions would allow m = 1 to dominate up to 
the point where strong non-linearities lead us to stop the 
simulations. 

In agreement with the results of figure El the modes we 
obtain have their corotation at « 1.1 — 1.2 tmso- Com- 
paring with the radial profiles shown in figure |B| we see 
that this does not correspond to an extremum of , but 
rather to a region where its slope is positive and large: 
this means that, rather than the RWI obtained with the 
free-flow boundary in section 13,31 we arc dealing here 
with an instability which is more of the AEI type (which 
relies on a reflecting boundary rather than on modes lo- 
calized in the extremum of although the effect of 
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the relativistic rotation curve forbids a clear separation 
between the two mechanisms. The strong slope of Lb, 
due to the relativistic rotation curve, explains the large 
growth rates we obtain. 

4. DISCUSSION 

In this paper wc have put together the physics of disko- 
seismologic models, of the AEI, and of the RWI, to find 
that normal modes are strongly unstable in the accretion 
disk of black holes, especially when the disk magnetiza- 
tion is of the order of equipartition. The instability mech- 
anism, relying on a coupling between spiral and Rossby 
waves in the disk, is the same that makes unstable the 
AEI and RWI, with variations depending on the condi- 
tions at the boundary between the disk and the black 
hole magnetospliere. 

We find that, if the gas is allowed to flow freely across 
the Marginally Stable Orbit, the dominant mode is a one- 
armed (to = 1) feature, not suitable to explain the high- 
frequency Quasi-Periodic Oscillations of microquasars. 
On the other hand wc have argued that, if the disk carries 
a net poloidal flux which is advected toward the central 
region, the complex magnctosphcric structure involved 
in the Penrose and Blandford-Znajek processes might 
act as a reflecting boundary for waves propagating in 
the disk. Using this assumption we find, from linear 
analysis and numerical simulations, that all the modes 
with TO > 2 have very similar pattern speeds {uj /m) and 
growth rates, while the to = 1 is less unstable. Small 
differences depending on additional physics {e.g. three- 
dimensional effects, a realistic equation of state, viscos- 
ity from smaller-scale MRI turbulence) can select any of 
them, which would thus appear very near multiples of an 
unseen fundamental frequency, as the HF-QPO do. This 
would readily explain the 2:3 (and sometimes higher) fre- 
quency ratio observed between these QPO. 
The modes are already unstable in non-magnetized 
disks, in which case they become forms of the RWI or 
Papaloizou-Pringle instabilities affected by the relativis- 
tic rotation curve. However they are much more unstable 
when the disk is threaded by a magnetic field approach- 
ing equipartition with the gas pressure (/3 > 1). The HF- 
QPO are observed in the Very High, or Steep-Power Law, 
state of the sources (see e.g. iMcClintock fc R.emillardI 
20061). nea r the jet line in the color-color diagrams of 
Fender et a l. (2004), except maybe in GRS 1915-1-105. 
This state is characterized by both a strong thermal 
disk emission and a steep power-law tail, probably from 
Comptonization in the corona. This contrasts with the 
Low-Frequency QPO, seen in the low-hard state where 
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the disk is hardly visible and the coronal emission is dom- 
inant, and that we consider as a manifestation of the AEI 
at /3 < 1. These modes share the property that the ac- 
cretion energy they extract from the disk is stored in a 
Rossby vortex, and we have shown ijVarniere fc Tagged 
l2002t) that a substantial fraction of this energy can be 
emitted toward the corona as an Alfven wave. Energiz- 
ing the corona by dissipation of this wave would thus 
provide a natural explanation for the spectral character- 
istics of the sources when these QPO are observed. 
Given these properties, it is natural to assume that the 
HF-QPO are seen at times where the magnetic field in 
the disk is below equipartition, allowing instability both 
at small scales by the MRI (explaining the thermal disk 
emission) and at large scale, in the vicinity of the MSO, 
by our instability mechanism. Further work should con- 
sider the relevance of this interpretatio n to the 'Magnetic 
Floods' scenario ()Tagger et al.l 1200^ we have recently 
presented for some of the cycles of GRS 1915-1-105. In 
this respect it is already remarkable that in this scenario, 
just as in the present work but from entirely different ar- 
guments, we were led to highlight the importance of the 
magnetic flux trapped between the inner disk edge and 
the event horizon of the black hole. 
On the other hand, this work should lead to take with 
caution estimates of the black hole spin, derived from a 
tentative identification of the QPO frequency with the or- 
bital frequency at the MSO, for a Kerr black hole whose 
mass is known. Here we find that the QPO frequen- 
cies are 2, 3, ...times a fundamental that is typically 
« .8 times the MSO rotation frequency. This factor 
may change in fully relativistic computations, or in ones 
where additional physics such as three-dimensional ef- 
fects, a realistic equation of state, viscous transport from 
small-scale turbulence, or the interaction with the central 
magnetospheric structure, would be taken into account, 
affecting in particular the radial profiles of density, tem- 
perature and magnetic field in the vicinity of the MSO. 
Much more elaborate simulations will be necessary to 
reliably assess these effects. 
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